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Abstract 



^L | The C 12 {C 12 ,a}Ne 20 fusion reaction rate is shown to be extremely sensitive to 

variations in the strong coupling constant, a s . Connection with a s is established by 
using the one boson exchange potential for the nuclear forces and then relating the 
meson masses and coupling constants of this model to QCD. The implications of a 
cosmological evolution of these constants on carbon burning in Type la Supernovae 
are discussed. 



1 Introduction 



Type la Supernovae (henceforth SNela) are thought to result from the ignition of a degen- 
erate Carbon/Oxygen stellar core. This occurs when the mass of the core increases until 
the energy released via C 12 +C 12 fusion can no longer be radiated away by neutrinos. SNela 
explosions are badly understood and it is not clear (for example) at what mass the core 
begins runaway burning, where in the core burning is instigated and what kinds of progen- 
itor systems lead to SNela's [0, §]. At low redshifts there is a strong relationship between 
SNela luminosity rise times and peak luminosities, rendering these objects potentially useful 
standard candles for cosmology. 

The low apparent luminosity of high redshift (0.3 < z < 1.0) supernovae can be interpreted 
as evidence for a FRW universe with a cosmological constant of A « 0.7 or as SNela 
evolution. It has been suggested that the peak luminosity of a SNela is higher if there 
is a lower Carbon/Oxygen ratio in the core [Q]. However, higher mass stars which evolve 
quickly and should therefore dominate the SNelA population at high redshifts have lower 
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Carbon abundance cores, creating a problem for this evolution argument. At the same time, 
recent studies suggest that the rise times of high redshift SNela's are shorter than the rise 
times of local SNela objects, indicating a possible evolutionary effect ||. The astrophysical 
analysis assumes that the nature of fundamental physical interactions have remained the 
same throughout the history of the Universe which, until tested rigorously, may not be the 
case. 

Over the past 20 years it has become clear that theories involving the presence of higher 
dimensions can provide a framework in which the four forces of nature may be described by 
a unified theory @. The fine structure constants in these theories are related to the size 
of the higher dimensions so it is neccesary to provide a mechanism to stabilise these sizes. 
As it is still not entirely clear what this mechanism may be, any observational evidence for 
the evolution of fine structure constants over cosmological time scales might shed light on 
the problem. At the same time, changing fine structure constants would lead to changes in 
the behaviour of supernovae. Most physical processes would be affected including nuclear 
fusion cross sections and energy yields so the evolution and luminosity of the explosion would 
certainly be different. In || the effect of the time variability of the gravitational constant on 
the Chandrasekhar mass and the background cosmology were calculated and used to reconcile 
the supernova observations with a A = open universe. In this article we investigate the 
effect of changing electromagnetic and strong fine structure constants on the carbon fusion 
responsible for SNela explosions. As i) it is not clear as to whether SNela explosions occur 
at the Chandrasekhar mass |2| and ii) it is not clear how sensitive explosion luminosities 
are to the mass of the progenitor ||, we will be less ambitious than the authors of || and 
simply calculate the effect of changing a em and a s on the carbon fusion cross section and 
energy yields (Q values). 

1.1 Existing Limits on Variations in a em and a s 

Webb et al. [||] placed constraints on the fractional variation of a em as being about 1 part in 
10 5 at a redshift of z=l. The variation of a s is more subtle as, because of well known QCD 
renormalisation effects on the strong coupling constant, a s becomes considerably larger at 
low energies. The 3 colour charge one-loop renormalization group equation for the QCD 
coupling constant is given by 

21 2 * + Uln(^\ (1) 
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where p has dimension of inverse mass or length and p c corresponds to some energy scale at 
which a s is known. This energy scale has to be within the regime of the three light quarks, 
up down and strange, in order that the above relation holds (at higher energies, virtual 
pairs of charm, bottom and top quarks will lead to corrections to the renormalisation group 
equation). To say what a em means today is straight forward, but the growing value of a s 
at low energies means that any QCD physics observed today occurs within highly bound 



systems such as the proton. Also, as one can see by examining the renormalization equation, 
small changes in a s at high energies correspond to very large changes in a s at low energies. 
Consequently, not only is it necessary to decide how much one is permitted to change a s 
but also at what energy scale this constraint is imposed. A very strict limit on the varying 
of a s has been provided by considering capture cross sections for slow neutrons [^Oj] . If one 



saturates this limit and assumes the rate of change is constant, a 2% limit on the fractional 
variation of a s is obtained for a time 10 10 years ago. Although these capture cross sections 
were observed for thermal energies much lower than the mass of the proton, the absorption 
will depend on the wave function of the absorbed neutron within the nucleus. This wave 
function will involve a potential derived from the interactions of nuclei at an energy scale of 
the order of the nuclear size and hence the proton mass. This is the energy scale at which 
a s has been varied in this study. 

2 The C 12 +C 12 Fusion Reaction 

2.1 Fusion Cross Section Theory 

The reaction rate for the fusion of two nuclei is given by p3 

/ 8 \ 1/2 1 r°° ( E \ 

where /i is the reduced mass, T the temperature, E is energy, a(E) is the reaction cross 
section, k is Boltzmann's constant and v is the velocity of particles in the medium. The 
cross section can be written 

a(E) = ^exp(-2n^^)S(E) (3) 

where Z\ is the number of protons in the target nucleus and Z 2 refers to the number of 
protons in the incident nucleus. S(E) parametrises the resonant structure due to different 
channels being favoured at different energies. Theoretical estimates of S(E) for C 12 +C 12 
fusion have been particularly controversial as early laboratory experiments suggested the 
presence of pronounced resonances at astrophysical energies (up to a few hundred keV) [ |TTf . 
These were thought to be due to the phenomenon of absorption under the barrier where the 
two nuclei form an intermediate nuclear molecular state. Several approaches to modelling 
this interaction were investigated, none of which predicted the observed resonant structure 
entirely successfully Jl2]| ]13[ Jl4| Jy| . However, in f36fl , it was shown that carbon fusion was 



modelled acceptably by the WKB approximation where the cross section is written \\t 



a(E) = n\ 2 (E)Y,(2l+m(E) (4) 



A being the reduced deBroglie wavelength of the incident nucleus and I being the angular 
momentum number. The transmission coefficients Ti(E) are written 

T l (E) = [l + exp(S l )}- 1 (5) 

where Si is the classical action evaluated between the inner and outer classical turning points, 
Rj and R Q . 
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S l = 2^J R °[V l (r)-E]dr (6) 

Here, [i is the reduced mass and the potential V/(r) is given by 

V5(r) = V coul (r) + V nuc (r) + ^ + *f ■ (7) 

The Gamow Peak energy, Ec a m, is the favoured energy for fusion in thermalised systems and 
corresponds to the energy at which the product of the Maxwell-Boltzmann distribution and 
the barrier penetration probability grows highest. Evaluation of the dominant £=0 channel 
cross section at Ea am will give a reasonable estimate of the change in the astrophysical fusion 
rate as the nuclear and Coulombic potentials change. The Gamow peak energy is given by 
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2.2 Inter-Nucleus Potential Model 

In order to evaluate the cross section, the potential due to both Coulombic and nuclear forces 
between the two carbon nuclei must first be calculated. This was done by modelling each 
carbon nucleus as a solid sphere of constant density with radius 

R = Voton^ 1/3 = 2.28984r proton . (9) 

Experimental data shows that this saturated density ansatz is a good approximation, al- 
though there are models thought to be more realistic where density tails off exponentially 



from the surface fl4j| . Equation (9) shows that one can calculate how the radius of any 
nucleus changes when constants are altered given that we know how the proton radius re- 
acts. The success of the saturated density model also indicates that the nucleus structure 
is dominated by short range effects and consequently by the strong force as opposed to 
electromagnetism. 

Next, we assume that the twelve units of nuclear charge (Baryon Number) and six units 
of Coulombic charge are evenly distributed throughout the nucleus. Then it is possible to 
do a numerical integration by approximating the two nuclei as spherical sections of a three 
dimensional grid. The potential energy at each individual site due to the nuclear charge 



located at one of the sites on the other nucleus is given by the single nucleon-nucleon one 
boson exchange potential with appropriate weighting. Because spins of individual nuclei are 
neglected in this simple model, the nucleon-nucleon potential is vastly simplified as we shall 
see in the next section. 

The cross section was calculated with all parameters fixed as observed in nature, then a em 
and a s were altered independently. When a em is altered, the Coulombic component of the 
inter-nuclear potential and the binding energies of the nuclei are altered. When a s is altered 
at some energy, the energy scale at which quark confinement occurs (approximately when 
a s (p) in (1) goes to one) also changes. This tells us how the proton radius and mass, and 
therefore the nuclear radius and mass, change and the change in the masses of the mesons in 
the boson exchange potential also can be calculated as we shall see later. Again, the binding 
energies and consequently the energy liberated by the reaction are changed by varying a s . 
These effects were included in the simulations. 

3 The One Boson Exchange Potential 

In order to calculate how the nuclear force between two C 12 nuclei will vary with QCD param- 
eters it is first necessary to consider the single nucleon-nucleon potential. This is described 
quite successfully by the one boson exchange potential (OBEP) P5j] . The most significant 
interactions of the four meson fields that play a role in the OBEP are the pseudoscalar in- 
teractions of the isovector tt, the attractive exchange of the scalar boson a, the repulsive 
vector interactions of the isoscalar u and the tensor interactions of the isovector p |35| . The 
tensor potential has an average value of zero upon integration over all angles of incidence 
| 37j so the p potential will be neglected as we will not be using a complicated C 12 model. 
Further simplification results from the fact that the pseudoscalar potential describing single 
7r exchange contains a factor o\ ■ a 2 taking into account the relative spins of the two spin 
1/2 nucleons. The value of this product is either -3 for the singlet (degeneracy 1) spin state 



or +1 for the triplet (degeneracy 3) spin states |38| so if we make no assumptions about the 
relative orientation of the spins of the nuclei in the model, the average force due to single 
Ti exchange vanishes. We are subsequently left with the scalar a and the repulsive vector uj 
potentials which can now be written as [[35] 
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Here, m u , m a are the masses of the omega and sigma mesons, g a NN, QujNN are the respective 
meson-nucleon-nucleon coupling constants, f u is the omega meson decay constant and m n is 



the nucleon mass. The next step is to relate the parameters of this potential to the strong 
fine structure constant. 

4 Meson Masses and Coupling Constants from QCD 

As we shall see later in this section, the QCD vacuum is characterised by quark and gluon 
condensates. In order to make contact between the OBEP and a s it is necessary to derive 
the meson masses and coupling constants in terms of these condensates. 

4.1 a Meson Mass and Coupling Constant 

QCD sum rules were used by Ioffe to obtain an approximate expression for the nucleon mass 
in terms of the quark condensate [IS| 

= -2(2,)^) 

M 2 v ; 

where m n is the mass of the nucleon and M is the Borel mass. This is a characteristic 
mass scale of the system introduced to regularise the series of contributions from higher 
dimensional operators which enter into light quark calculations |27j • By assuming the Borel 
mass is the nucleon mass, one can write the approximate equation 

m n = [-2(2n) 2 (qq)} 1 / 3 (13) 

However, nucleon mass can be explained more satisfactorily by the Nambu Jona-Lasinio 



(NJL) model [19| [20] . In this model, the nucleon acquires mass by interacting with the 
quark condensates of the vacuum via a scalar field. Brown and Adami use an argument 
based upon the dependence of both sides of Ioffe' s expression on N c , the number of colour 
charges in the theory (3 for the SU(3) of QCD) to identify the scalar field of the NJL model 
with the sigma meson. This results in the identification of the Borel mass not with the mass 



of the nucleon but with the sigma mass pl|, [p^| . Using a constituent quark model for this 
process, they were able to write 

m n = -|%<g<z). (14) 

3 m A a 

Comparison of the above equations completely defines the sigma meson mass in terms of the 
quark condensate and the nucleon mass 

m c = g aNN (15) 

\3 m n J 

and sets the coupling constant value 

g 2 aNN = 12vr 2 . (16) 

Consequently, the change in the proton radius (inverse mass) will have to be taken into 
account when calculating the change in the a meson mass. 
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4.2 uj Meson Mass and Coupling Constant 

The mass of the u meson was predicted very accurately as long ago as 1977 when Chan 



considered the quark- loop expansion of the vector meson form factor |23 ]. The result is 

simply 

3.1/2 



m w = m p 



2^/,(-) . (17) 



The degeneracy between the u and p meson masses is broken by the fact that each has 
different total isospin and therefore mixes with virtual meson pairs differently (see e.g. ||24|| ). 
However the approximation of Chan is still in excellent agreement with observation (about 
two percent [p^1,[p5|). Now the omega coupling constant is related to the omega mass and 



the pion decay constant, f w by |25 



fc ™=S0f=*^ (i8) 

determining all the parameters of the model in terms of the pion decay constant and the 
quark condensate. 

4.3 The QCD Vacuum 

A classical solution of the equation of motion in a Euclidean space-time sometimes represents 



quantum tunnelling through a barrier in a Minkowski space-time [|PT |. Consequently, the 
discovery in the 1970's of an exact finite action Euclidean solution to the classical Yang- 
Mills equations lead to the realization that the QCD vacuum is highly non-trivial. At short 
distances where the effective coupling (oc tJoT s ) is small, the dominant contributions to the 
path integral characterising the true vacuum are from tunnelling transitions between an 



infinite number of degenerate vacua [pOfl . The vacuum is therefore determined at these 
distances by instanton effects rather than perturbation theory. The way the density of 
instantons changes at different length scales p in the colour SU(3) theory is given by the 
weighting function [BO 



Dsmip) = (0.10) (JL) e X p (-JL) . (19) 

As one goes to larger distances, the effective coupling rises and so one must consider quantum 
corrections to the vacuum. However, even though individual tunneling amplitudes are small, 
the rapidly increasing number of instantons at larger radii means that tunneling contributions 
are still more important than loop corrections. The qualitative explanation for confinement 
in this model is that small instantons and anti-instantons act to screen larger instantons, 
the effect of which are only felt at distances corresponding to their size, and the effective 
coupling rises rapidly at larger distances - hence confining any particle with colour charge. 



Quarks carry colour charge so the distance at which the coupling constant becomes large is 



directly associated with the proton radius [32 



The fact that the classical description forms the bulk of the vacuum structure at small 
radii allows the use of the dilute instanton gas approximation to determine the self interaction 



of gluons and consequently an expression for the gluon condensate [p7 
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where G^ u is the gluon tensor and p c is an infrared cut-off corresponding to the length scale 
at which the effective coupling goes to one. 

In the above picture, the QCD vacuum energy is obtained from the trace of the stress 
energy tensor which receives contributions from the quark and gluon condensates and any 
non-zero bare quark masses. The QCD vacuum structure can also be obtained by considering 
the effective potential for composite operators j25fl . In this scheme, the vacuum energy will 
be a function of the bare quark masses and the quark condensate alone. By comparing these 
two methods it is possible to relate the gluon and quark condensates. In this work, we do 
calculations at the chiral limit where the up and down quarks have no bare masses. Quark 
masses are hence entirely dynamically generated by vacuum screening effects. At this limit, 



the relationship between the quark and gluon tensor has been recently established J29] 
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where e 1 can be shown to be |E9 
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and (3{a s ) is the perturbative 2 loop beta function 
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N c and nj being number of colours and flavours respectively. 
Equation (20) can be written as an expansion in a s (p c ) J32"H 
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the dominant contribution to the expansion coming from the unity term inside the square 
brackets. The value of the gluon condensate can be fixed by observation of the mass splitting 
of the ground states of J/ip and r\ c mesons WA - By using the observed value of 0.012 GeV 4 



and setting a s (Mp rot ) = 1 one can obtain a value of p c ~ l/125MeV, within a factor of 8 of 
the proton radius. We then recycle this factor 8 to replace p c with the inverse proton mass 
in equation (24). As we know how the proton mass varies with the strength of the strong 
force, we assume that p c varies in the same way. 

The dynamically generated quark mass in the chiral limit at low momentum scales is 
given by |29| 

m dyn = (-qq) 1/3 (25) 

then the pion decay constant can be written |26| 

U = 7f^I7I m ^ = 0.34025(-gg) 1 / 3 (26) 



which determines all of the meson parameters in terms of quantities from QCD. 

5 Binding Energies and Q- Values 

Here we make the assumption that despite the fact that shell closure effects render the liquid 



drop model inaccurate for nucleons with A<20 pBfl , the model can give an indication of the 
dependance of the nuclear binding energy on a em and a s . The Semiemperical Mass formula 
is given by 

Z 2 (Z — N) 2 
M(A) = ZM prot + NM neut - ai A + a 2 A 2/3 + a 3 —j, + a 4 - —^- + 5(A) (27) 

where the terms with coefficients oi - a^ are the volume, surface, Coulomb and asym- 
metry terms respectively. The dominant channel for Carbon Burning is C 12 {C 12 ,a}Ne 20 
(C 12 +C 12 =^Ne 20 + a). The volume contribution to the Carbon-Carbon system is the same 
as that for the a-Neon system and consequently cancels out of the Q- value calculation. The 
pairing term 5(A) gives a positive contribution to the mass for nuclei with odd numbers of 
both protons and neutrons (odd-odd), a negative contribution of the same magnitude for 
even-even nuclei and zero contribution for odd-even nuclei. As Carbon nuclei, a particles 
and Neon nuclei are all even-even, the final state has the same ammount of pairing energy 
as the inital state and we do not need to calculate the dependance on r prot of the binding 
energies for this reaction. The assymetry term can be neglected for all three nuclei, leaving 
us with the surface, Coulomb and pairing terms. At this stage we still label the strength of 
electromagnetism directly by a em but the contrived root from a s (p) to observable physics 
makes it easier to signify the strength of the strong force by r prot . The empirical prefactors 
for the surface and Coulomb terms can be estimated quite well by using a simple Fermi gas 
model for the nucleus. 

This model uses the fact that only two neutrons and two protons can be in the same 
place with the same momentum to relate the density of nucleons to the highest occupied 



momentum state, the Fermi momentum kp. 
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Since kp is proportional to the constant nuclear density and the nucleon mass is the inverse 



of the nuclear radius, kp ocr™./ t . The Fermi energy is therefore given by ep ocr pro ' t . In [R9 
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it is shown that the surface term coefficient as calculated by Fermi Gas considerations is 
only about 3% different from the empirical value. The surface coefficient is obtained from 
the model by remembering to neglect any unphysical momentum states which occur when 
one or more components of the momentum vector are zero. The result is 

9vrep _ 4/3 

«2 = -T7T, ocr U . (29) 

4:0k F r prot 

Similarly, the Coulombic part of the energy can be derived from the Fermi gas model, 
although not as accurately as the surface term. The full expression for the normal and 
altered terms is then 
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The Q value for the C 12 {C 12 ,a}Ne 20 reaction using the semi empirical mass formula is 
hence 

Q = 0.595a 2 ( SU rface) ~ 7 ■^^ a 3(Coulomb) (31) 

This yields a value of 5.37 MeV, about 15% out from the observed value of 4.62 MeV |34]] but 
we were expecting a discrepancy, as the liquid drop model is not entirely reliable for A<20. 
We shall continue using this approximate method in the hope it will give us an indication of 
the magnitude of the effect of changing fine structure constants. 

6 Results and Discussion 

When a em is varied in the equations set out in section 2, the Carbon fusion cross section 
changes significantly more. Saturating the bounds set by Webb of lxlO -5 at a redshift of 
z=l, it is found that the fractional cross section change is 5.4xl0 -4 , a factor of 50 greater 
fractional change than a em . However this is still very small and will make probably make 
very little difference to overall luminosities as we shall see below. 

Using the equations in section 4 it was found that for a 2% increase/decrease in a s at the 
proton mass, the meson masses increase/decrease by about 2.5% (both are oc (qq) 1 ^ 3 )- The 
same change in a s resulted in an approximate 1.2% change in the radius of the proton, the 
radius decreasing as the coupling was increased. Next, the effect of the changing potentials 
on the the 1=0 cross section at the Gamow peak corresponding to a temperature of lx 10 9 K 
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were calculated (lxl0 9 K is the typical temperature quoted for the time when the supernova 
is making the transition from quiescent to explosive burning [J40"f). 

Increasing a s (M Prot ) by 2% was found to decrease the cross section by 51% and decreasing 
a s (Mp rot ) by 2% increases the cross section by 115%, so both ways change the cross section 
by about a factor of two. The Q-value of the reaction also rises as a s goes down, with 
AQ/Q « 0.04. 

The fact that a 2% change in a s yields can halve or double the fusion cross section is 
remarkable and could lead to interesting effects. As an example, one might expect that if 
a s had been 2% smaller at some time since the big bang, Carbon burning would take place 
much more readily, meaning that degenerate cores would explode at lower temperatures - 
hence lower masses. A lower mass core would probably result in a lower luminosity explosion 
iTi| ] and we would have a new interpretation of the apparent low luminosity of high redshift 



SNela's. Unfortunately, things are not this simple because of the extreme temperature 
sensitivity of the Carbon fusion reaction. 

The apparently large change in the fusion cross section would be compensated for by a 
very small change in temperature. One way of quantifying this effect is by taking the fusion 
cross section and multiplying it by the Q value to obtain a quantity signifying the energy 
yield for the three different values of a s . It turns out that a 6.7% increase in temperature 
would compensate for a 2% decrease in a s . In a SNela this would certainly result in a change 
in the rate at which the transition from quiescent to explosive burning occurs and the peak 
luminosity would also change. However, the exact way this would effect the relationship 
between rise-time and peak luminosity is difficult to quantify without a more complete 
understanding of SNela's. 

Models of Type la Supernovae are extremely complex, using detailed nuclear reaction 
schemes that take into account many different cross sections. The fact that these models 
are still not entirely succesful after thirty years of endeavour show how complicated the 
mechanisms occuring in SNela's are. It seems there are many other effects which may be 
responsible for the apparent evolution of SNela luminosity which do not involve speculative 
physics. However, if there is reason to believe in a cosmological evolution of a s , SNela's may 
be a good place to start looking. 
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